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calculations at hadron colliders, we derive the full set of initial-initial double real radiation 
antenna functions in integrated form. 

Keywords: QCD, Jets, Colhder Physics, NLO and NNLO Calculations. 



Contents 



1. Introduction 1 

2. Initial-initial antenna subtraction at NNLO 3 

3. Antenna functions for double real radiation 5 

3.1 Definition of the antenna functions 5 

3.2 Phase space factorisation and mappings 7 

4. Integration over the double real radiation phase space 9 

5. Results 11 

5.1 Quark-antiquark antenna functions 11 

5.2 Quark-gluon antenna functions 13 

5.3 Gluon-gluon antenna functions 15 

6. Conclusions 16 

A. Master integrals 17 

A.l Hard region 19 

A. 2 Collinear xi region 25 

A. 3 Collinear X2 region 29 

A.4 Soft region 33 



1. Introduction 

Jet production at high energy colliders allows to reconstruct the parton-level dynamics of 
the hard interaction processes. Jet cross sections axe therefore ideally suited for precision 
studies of QCD [1] to determine the strong coupling constant or parton distribution func- 
tions. They can be measured to a very high experimental accuracy, which has to be met 
by an equally high theoretical precision in order to perform meaningful comparisons. This 
precision demands the inclusion of ncxt-to-ncxt-to-leading order (NNLO) QCD corrections. 

NNLO calculations of observables with n jets in the final state require several ingre- 
dients: the two-loop corrected n-parton matrix elements, the one-loop corrected (n -|- 1)- 
parton matrix elements, and the tree-level (n -|- 2)-parton matrix elements. For many 
massless jet observables of phenomenological interest, these matrix elements are available 
for some time already. 

The (n+ l)-parton and (n + 2)-parton matrix elements contribute to n jet observables 
at NNLO if the extra partons are unresolved or are clustered to form an n-jet final state. 
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Consequently, these extra partons are unconstrained in the soft and collinear regions, and 
yield infrared divergences. In these cases, the infrared singular parts of the matrix elements 
need to be extracted and integrated over the phase space appropriate to the unresolved 
configuration to make the infrared pole structure explicit. The single soft and collinear 
limits of one-loop matrix elements [2-6] and the double unresolved limits of tree-level matrix 
elements [7-10] are process-independent, and result in a factorisation into an unresolved 
factor times a matrix element of lower multiplicity. 

To determine the contribution to NNLO jet observables from these configurations, one 
has to find subtraction terms which coincide with the full matrix element in the unresolved 
limits and are still sufficiently simple to be integrated analytically in order to cancel their 
infrared pole structure with the virtual contributions. Often starting from systematic 
methods for subtraction at NLO [11-14], several NNLO subtraction methods have been 
proposed in the literature [15-22], and are worked out to a varying level of sophistication. 

For observables with partons only in the final state, an NNLO subtraction formal- 
ism, antenna subtraction, has been derived in [23]. The antenna subtraction formalism 
constructs the subtraction terms from antenna functions. Each antenna function encapsu- 
lates all singular limits due to the emission of one or two unresolved partons between two 
colour-connected hard radiator partons. This construction exploits the universal factori- 
sation of matrix elements and phase space in all unresolved limits. The antenna functions 
are derived systematically from physical matrix elements [24]. This formalism has been 
applied in the derivation of NNLO corrections to three-jet production in electron-positron 
annihilation [25,26] and related event shapes [27,28], which were used subsequently in 
precision determinations of the strong coupling constant [29-34]. The formalism can be 
extended to include massive fermions [35]. Antenna subtraction has also been used as the 
starting point for a parton shower algorithm [36], using the fact that antennae capture the 
unresolved limits of QCD amplitudes to achieve a leading-log resummation. 

For processes with initial-state partons, antenna subtraction has been fully worked out 
only to NLO so far [37]. In this case, one encounters two new types of antenna functions, 
initial-final antenna functions with one radiator parton in the initial state, and initial-initial 
antenna functions with both radiator partons in the initial state. The framework for the 
construction of NNLO antenna subtraction terms involving one or two partons in the initial 
state has been set up in [38,39] in the context of a proof-of-principle implementation of the 
contribution of the purely gluonic contributions to di-jet production at hadron colliders. 
The initial-final and initial-initial antenna functions appearing in the NNLO subtraction 
terms are obtained from crossing the final-final antennae. Their integration has to be 
performed over the appropriate phase space. In the case of the initial-final antennae, 
this has been accomplished in [40]. The initial- initial one- loop antenna functions were 
integrated in [41] . Partial results have been obtained for the integrated initial-initial double 
real radiation tree-level antenna functions in [42] , where all antennae involving a secondary 
quark-antiquark pair were integrated. It is the aim of the present paper to complete 
the NNLO antenna subtraction scheme for hadron collider processes by integrating the 
remaining initial-initial double real radiation antenna functions. 

Other approaches to perform NNLO calculations of exclusive observables with initial 
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state partons are the use of sector decomposition (partly in combination with subtraction) 
and a subtraction method based on the transverse momentum structure of the final state. 
The sector decomposition algorithm [43] analytically decomposes both phase space and loop 
integrals into their Laurent expansion in dimensional regularisation, and performs a sub- 
sequent numerical computation of the coefficients of this expansion. Using this formalism, 
NNLO results were obtained for Higgs production [44] and vector boson production [45] 
at hadron colliders. Both reactions were equally computed independently [46] using an 
NNLO subtraction formalism exploiting the specific transverse momentum structure of 
these observables [20] , which was also applied to compute NNLO corrections to associated 
WH production [47] and photon pair production [48] . A very promising approach could be 
the combination of subtraction with sector decomposition [21], which was recently applied 
in the calculation of NNLO corrections to the quark-induced processes in top quark pair 
production [49]. 

This paper is structured as follows: in Section 2, we briefly summarise the application 
of the antenna subtraction formalism to hadronic collisions, which was developed in [37- 
39,42]. The antenna functions required for double real radiation and their associated phase 

space mappings are described in Section 3 while the phase space integration is described 
in Section 4. The resulting integrated initial-initial antenna functions being too lengthy to 
be presented here, only their leading pole parts are given in Section 5. The full results are 
attached separately to this paper as a FORM [50] file. Section 6 contains our conclusions 
and an outlook. The newly derived phase space master integrals required for our calculation 
are collected in the Appendix. 

2. Initial-initial antenna subtraction at NNLO 

The partonic double real contribution to an NNLO m-jet cross section reads 

d^NNLO = ^NNLO d$„+2 {ki, . . . ,km+2]Pl,P2) 

perms '^"^+2 

X |Mm+2(pi,P2;A;i,...,A;^+2)P J^+^HA;i,...,A;„+2;Pi,P2) • (2.1) 

In this equation, \M.m+2{pi-,P2\ki^ ... , km+2)\^ stands for the colour-ordered 2 ^ m + 2 
amplitude norm squared (or, for subleading colour contributions, for the appropriate sum 
of interference terms between colour-ordered amplitudes). The symmetry factor 5*^+2 ac- 
counts for identical partons in the final state and ^p^j-j^g denotes the sum over all configu- 
rations with m + 2 partons. The next-to-next-to leading order normalisation factor M§^j^q 
includes all QCD-independent factors as well as the dependence on the renormalised QCD 
coupling constant a^. It is related to the normalisation factor present at leading order, 
Nlo, which depends on the specific process and parton channel under consideration by, 

^NNL0-^fL0[-^'j (2.2) 



where. 
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C(6) = (47r)^^, (2.3) 
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C{e) = (47r)^e-^^. (2.4) 

The initial state momenta are labelled as usual as pi and p2 whereas the m+2 momenta 
in the final state are labelled ki,... , km+2 ■ The 2 — >■ m + 2 particle phase space is denoted 
as 

d$„j+2(A;i, . . . , km+2;pi,P2) = [dki] . . . [dkm+2] (27r)'^5'^(pi +p2 - ki - . . . - km+2) (2.5) 

where we have introduced the abbreviation [dki] = d'^kiS'^ {kf) / (p.ny^^ . The jet function 
Jm ^ {kl , ■ ■ ■ , kfn+2 ;Pi,P2) ensures that out of (m + 2) final state partons, an observable 
with m jets is built. The incoming parton momenta pi,p2 serve as reference directions to 
define transverse momenta and rapidities of the jets. 

The real radiation contribution (2.1) contains infrared divergencies, which arise when 
one or two of the final state partons become unresolved. The numerical phase space 
integration in (2.1) can therefore not be carried out without prior introduction of a regulator 
or a subtraction. A subtraction term da^jy^Q is defined on the same phase space as 
d^N^LO^ and it approaches its integrand in all unresolved limits. Consequently da^^^^Q — 
d^NNLO finite and the integration can be performed numerically. d&^j^j^Q is integrated 
over those parts of the phase space that contain infrared singularities and then combined 
with the virtual NNLO and the mass factorisation counterterm contributions to achieve 
the cancellation of infrared divergencies. 

The double real radiation contribution da^^j^Q can become singular if either one or 
two final state partons are unresolved (soft or collinear). Consequently, when constructing 
the corresponding subtraction term daf^j^^Q in the antenna subtraction method, we must 
distinguish the following configurations according to the colour connection of the unresolved 
partons: 

(a) One unresolved parton but the experimental observable selects only m jets. 

(b) Two colour-connected unresolved partons (colour-connected). 

(c) Two unresolved partons that are not colour connected but share a common radiator 
(almost colour-unconnected) . 

(d) Two unresolved partons that are well separated from each other in the colour chain 
(colour-unconnected) . 

(e) Large angle soft gluon radiation. 

All cases except (b) and (e) can be accounted for by combinations of three-parton 
antenna functions, which were derived and integrated already in the context of antenna 
subtraction at NLO. Case (e) can always be described [39] by soft antenna functions in 
final-final kinematics, that were derived in [25, 26]. In case (b), tree-level four-parton 
antennae are needed. The corresponding subtraction term for both radiators in the initial 
state reads: 

d^Nirio = ^NNLO d^rn+2{ku km+2; Pi , P2) 

m+2 '^"^+2 
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X XI i^hk - ^l,jkXiL,K - ^i,kAj) \-^m{XiPi, Xipi; ki,..., k, k,---, W2)P 

il jk 

X Jin'\h,---,ki,ki,...,k^j^2;xiPi,xipi) . (2.6) 

The sum runs over all colour-adjacent pairs j, k and implies that the hard momenta i, I are 
chosen accordingly. Tree-level antenna functions are generically denoted X^, we will give 
the precise definition in Section 2. Xfi -^^ contains all the unresolved limits associated with j 
and/or k becoming becoming unresolved between the initial-state emitters i and I. Its single 
unresolved limits are subtracted by products of three-particle antenna functions, ensuring 
that dcT^^'^Q is only active in the double unresolved region. The reduced matrix element 
Ai is evaluated with a set of m on-shell momenta which are obtained from the original 
ones by rcscaling the initial-state momenta pi and pi, omitting the final-state momenta kj 
and kk and applying a Lorentz boost to the remaining final-state momenta [37]. 

To analytically integrate the subtraction term (2.6), we employ dimensional regulari- 
sation in d = 4 — 2e dimensions. We use the factorisation of the phase space 

d^m+2(fel, • • • , km+2\Pl,P2) = d$m(^l, ■ ■ ■ ,ki,ku ■ ■ ■ , ^m+2; XiPi,X2P2) 

X(5(xi - xi) 5{x2 - X2) [dfej] [dkk] da;ida;2 , (2.7) 

where xi and X2 are defined throughout the phase space. In the case of coUincar initial 
state radiation, they are identified with the coUinear momentum fractions of the partons 
entering the hard scattering process. 

The kinematics of the reduced matrix element appearing in (2.6) depends on the 
mapped momenta, and thus consequently on q^,x\,X2- The integration of the antenna 
functions must therefore retain the dependence on xi,X2- The integrated initial-initial 
antenna functions are defined as: 

^il,jk{xh XI, e) = l^^p J [dfefc] Xi XI Sixi - Xi) 5{xi - xi) Xf^^^ (2-8) 
where C(e) is given in eq. (2.3). 

3. Antenna functions for double real radiation 
3.1 Definition of the antenna functions 

Antenna functions are characterised by their parton content and their radiators, i.e. the 
pair of hard partons which they collapse to in the unresolved limits. Accordingly, we group 
them into quark-antiquark, quark-gluon and gluon-gluon antennae. They are derived from 
physical matrix elements associated to the decay of a colourless particle into partons [24]. 
The trcc-lcvcl antenna functions are obtained by normalising the three- and four-parton 
tree-level colour sub-amplitudes squared to that of the basic two-parton process: The final- 
final three- and four-particle antennae are respectively defined by: 

^ijk — '^ijk,IK |_yV^0^|2 ' 
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quark-antiquark antennae 

quark-gluon antennae 

^3° Dl{q,g,g),Dl{q,g,g) 
El El{q,q',q'),El{q,q',iy* 

gluon-gluon antennae 

F^{g,g,g) 
G's Gl{g,q,%,G'i{g,q,^y* 



Table 1: The tree-level three-particle initial-initial antennae and their distinct crossings. is 
symmetric under the interchange of the two gluons, Fg under cyclic interchange of its arguments. 

I A/fO |2 

^ijkl - ^ijkl,IL |^o^|2 ■ K'^-^) 

where S denotes the symmetry factor associated with the antenna, which accounts both 
for potential identical particle symmetries and for the presence of more than one antenna 
in the basic two-parton process. It is chosen such that the antenna function reproduces the 
unresolved limits of a matrix element with identified particles. The initial-initial tree-level 
three- and four-parton antennae denoted by X?^ ^ and ^.^ are obtained by crossing two 
partons to the initial state, starting from the final-final antennae. The distinct crossings 
of the three-parton antenna functions are listed in Table 1. Crossings which are free of 
unresolved limits arc marked with a **. 

The four-particle trcc-lcvcl antenna functions arc not determined by the species of the 
particles alone but also by the colour-connection. We distinguish leading-colour antennae, 
denoted by letters without tilde, where the particles are colour-connected in the order they 
are listed and subleading colour antennae, denoted by letters with tilde, where the gluons 
are photon-like. This notation has been established in [23,40]. The unresolved limits of the 
initial-initial antennae can be obtained from those of the final-final antennae by crossing. 
The crossing of the triple-collinear splitting functions is explained in [51]. 

Any two particles of a four-particle final-final antenna can be crossed to the initial 
state to obtain an initial-initial antenna; therefore one final-final four-particle antenna gives 
rise to six initial-initial antennae. Due to symmetries, at most four of these initial-initial 
antennae are different. The independent crossings [42] are listed in Tables 2-4. To make the 
colour connection clear, in this list we write out the arguments of the antennae explicitly, 
i.e we write X4 (l^,j,k,l^ (where T denotes an incoming particle) instead of Xf^^-^. The 
initial-initial antennae which are free of singular limits are not needed for the construction 
of subtraction terms, but their integrated form could prove useful for cross-checks. 
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quark-antiquark antennae 



^4 Al{q,g,g,q), {q, g, g, q) , {q, g,g,q), Al{q, g, g, q) 

A4 Al{q,g,g,^), Al {q, g, g, q) , Al{q, g, g, q) 

Bl B2{q,q',q',^, B2{q,q',q',q), Bl{q,q' ^qf* 

C2 C2 [q, q,q,'g),C2 [q, q, q, q) , (g, ?, q, q) ** , Cl{q,%q,^) 



Table 2: The tree-level four-particle quark-antiquark antennae and their distinct crossings. is 

symmetric under the interchange of the two photon-like ghions. The nonidentical-flavour antenna 
i?4 is separately symmetric under interchange of q' with q' and of q with q. 



quark-gluon antennae 







{q, 9,9,9), Dl {q, g,g,g), Dl {q, g, g, 


9) 


, Dl{q,g,g,g) 




El 


El\ 


[q, q', q', g) , El {q, q', q', g) , El (g, q>, 


1', 


g), El{q,q',q', 


9) 


El 


El\ 


[q, q', q', g) , El {q, q', q' , g) , El {q, q', 




g), El{q,q',q'., 


9) 



Table 3: The tree-level four-particle quark-gluon antennae and their distinct crossings. Due to the 
cyclic colour connection, is symmetric under interchange of the second and fourth gluon. 



gluon-gluon antennae 

F! F2{g,g,g,g),F2{g,g,g,g) 

G4 Gl {g, q,q,g),Gl {g, q, q, g) , Gl {g, q,%g),Gl {g, q, f , g) 
Gl Gl [g, q, q, g) , Gl [g, q, q, g) , Gl{g, q, f , g) 
Hi Hl{q,%q',q'),Hl{q,q,q',q') 



Table 4: The tree- level four-particle gluon-gluon antennae and their distinct crossings. F4 is 
symmetric under cyclic interchange of its arguments. G4 is symmetric under the interchange of 
the two photon-like gluons as well as under the interchange of q with q. has three symmetries, 
q ^q, q' ^q' and the flavour renaming q -n- q' , q -^r^ q' . 

3.2 Phase space factorisation and mappings 

The construction of subtraction terms requires a mapping from the original set of momenta 
onto a reduced set. The mapping interpolates between the different soft and collinear lim- 
its which the subtraction term regulates. An appropriate mapping for the initial-initial 
case, both for single and double unresolved configurations, has been discussed in [37]. By 
requiring momentum conservation and phase space factorisation, the phase space mapping 
is strongly constrained. The remapping of initial state momenta can only be a rescaling, 
since any transversal component would spoil the phase space factorisation. For two unre- 
solved partons j and k, a complete factorisation of the phase space into a convolution of 
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an m-particle phase space depending on redefined momenta only and the phase space of 
the unresolved partons j and k can be achieved with a Lorentz boost. This boost maps 
the momentum q = pi + P2 — kj — kk , with > and pi,P2 being the momenta of the 
hard emitters, into the momentum q = xipi + X2P2 , where xi and $2 are fixed in terms 
of the invariants as follows: 

.^^(su-s^2^l_y^ ^^^r s,2-s,,-s,, i_y 

\S12-Sij-Sik S12J \Sl2-Sj2-Sk2 S12/ 

These two definitions guarantee the overall momentum conservation in the mapped mo- 
menta and the correct soft and collinear behaviours . The two momentum fractions f 1 and 
X2 satisfy the following limits in double unresolved configurations: 

1. j and k soft: xi — )■ 1, X2 — > 1, 

2. j soft and k^ = zipi: xi ^ 1 — zi, X2 —> 1, 

3. kj = zipi and kk = Z2P2- xi ^ 1 - zi, ^2 ->■ 1 — Z2, 

4. kj = zipi, kk = Z2P1: xi ^ 1 - zi — Z2, X2 1, 

and all the limits obtained from the ones above by the exchange of p\ with p2 and of kj 
with k^. The construction of NNLO antenna subtraction terms also requires that all single 
unresolved limits of the four-parton initial-initial antenna functions X^j^, with radiators 
i and I, have to be subtracted, such that the resulting subtraction term is active only 
in its double unresolved limits. A systematic subtraction of these single unresolved limits 
through products of two three-parton antenna functions can be performed only if the NNLO 
phase space mapping turns into an NLO phase space mapping in its single unresolved 
limits. A detailed discussion of the corresponding translation between these two momentum 
mappings can be found in [37]. 

The factorisation of the (m -|- 2)-parton phase space into an m-parton phase space and 
an antenna phase space involving the unresolved partons j and k given in eq. (2.7) can 
equivalently be written as 

d^m+2{kl, . . . , km+2;Pl,P2) = d^mih, ■ ■ ■ , ^i, ki, . . . , Ai„i+2; XiPi,X2P2) 

X J 5{q^ - xiX2Si2)S{2{x2P2- xipi) ■ q) (3.3) 
X [dkj] [dkk] dxi dx2 , 

where J is the Jacobian factor defined by 

J = S12 (a;i(si2 - Sij - Sik) + X2(si2 - S2j - S2fc)) • 

Using this phase space parametrisation, we can express the integrated initial-initial 
antenna functions (2.8) as: 

Xljk{x,,X2,e) = J [dkj][dkk] J X, X2 d{Ci)S{C2)Xl^,, (3.4) 
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where 



Ci = - xi X2 si2 , 

C2 = 2{x2P2 - xips) ■ q . (3.5) 

In the following section, we describe how the integrals (3.4) are carried out by means of 
a reduction to phase space master integrals with constraints (3.5) and a derivation of the 
master integrals from their differential equations. 

4. Integration over the double real radiation phase space 

The initial-initial antenna functions have the scattering kinematics 

P1+P2 ^ kj + kk + q 

where q is the momentum of the outgoing colourless particle. The momenta satisfy: 

Pi=P2 = k] = kl = 0, q^ = f = xi X2 S12 • 

The four-parton initial-initial antennae need to be integrated over the phase space of the 
unresolved partons j and k. This integration yields a result which depends only on q^, x\ 
and X2- The dependence on is only multiplicative, according to the mass dimension of 

the integral. 

Using the optical theorem, these phase space integrals are expressed as cuts of forward 
scattering two-loop diagrams, where the J-function conditions fixing xi and X2 are intro- 
duced through non-standard propagators [52,53]. The full set of propagators appearing in 
the four-parton antenna functions is: 





= iPi - kjf 






Dki 


= ipi - kkf 






Dj2 


= {P2 - kjf 






Dk2 


= {P2 - kkf 








= {kj + kkf 






Djki 


= (Pi - kj - 


kkf 


! 


Djk2 


= (P2 - kj - 


kkf 






= (Pl +P2- 


k,f 




Dkl2 


= (Pl +P2- 


kkf 






= kj , 






Dk 








Djkl2 


= (Pl +P2~ 


kj — 


kkf - 


Djkl23 


= {P3 +P1 + 


P2 - 


kj ~ kk) 



where ps = X2P2 — xipi. The cut propagators are Dj, Dk, Djku, Djki23- Only integrals 
in which all cut propagators appear in the denominator are non- vanishing. Besides the 
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four cut propagators, the phase space integrals of the antenna functions contain at most 
four further propagators. Any scalar product involving the integration momenta can be 
written as a linear combination of an appropriately chosen set of seven linear independent 
propagators (including the four cut propagators). Using this, we are able to express all 
integrands appearing in the phase space integration of the antenna functions in terms of 
a linearly independent set of seven of the above propagators. For each of these sets, we 
derived relations between different integrals based on integration-by-parts identities [54], 
which arc solved by using the Laporta algorithm, using the implementation in [55]. As a 
result, all phase space integrals can be expressed in terms of 20 independent phase space 
master integrals. Of these, 10 were already computed in [42]. 

The set of master integrals which we denote by Ii{xi,X2, e) are functions of xi, X2 and 
e. We begin by factoring out the leading behaviour of the master integrals Ii{xi,X2,e) in 
the limits xi — ^ 1 and X2 1, keeping the exact e-dependence: 

Ii{xi,X2,e) = {l-xir^-^' {1-X2r'-^' J'i{xi,X2,e). (4-2) 

The integers mi, 1712 are characteristic to each master integral. The functions ^-^(xi, X2, e) 
are regular at xi = 1, at ^2 = 1, and at a:; 1 = 0:2 =1 and can be calculated as Laurent 
series with, at most, second order poles in e. 

The integrated antennae given by A'(xi, X2, e) are linear combinations of these master 
integrals Ii{xi, X2, e), with coefficients containing poles in e, as well as in (1 — xi) and 
(1 — X2). After the masters have been inserted into the integrated antennae, those take the 
form 

X{xi,X2,e) = (l-Xl)-^-2^(l-.^■2)-'-2^7^(xl,X2,e), (4.3) 

where 7?.(xi,X2,e) is regular at the boundaries xi = 1, X2 = 1, and at xi = X2 = 1. The 
e-expansion of the singular factors (1 — Xj)^^"^*^ is done in the form of distributions: 

(1 - x,)-i-2^ = -Ls{l-Xi) + J2 ^-^T^nixi) , (4.4) 

n 

with 

To evaluate the integrated antennae, we decompose the phase space into four regions 
depending on the values of Xi and X2. Those regions are given by: 

• xi 7^ 1, X2 7^ 1, which we refer to as the hard region 

• xi = l, X2 7^1, and xi 7^ 1, X2 = 1, referred to as collinear regions 

• xi = 1, X2 = 1, which we denote the soft region. 

In the hard region (xi 7^ 1, X2 7^ 1), harmonic polylogarithms of weight two appear in the 
0{€^) term of TZ. Therefore, the e-expansion of the master integrals in the hard region is 
needed at least up to the order at which terms of transcendentality two appear. 
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In the collinear regions (xi = 1 or X2 = 1), since the expansion in distributions in 
eq. (4.4) generates additional 1/e factors, the function TZ is required up to 0{e), where 
harmonic polylogarithms of weight 3 appear. The masters evaluated in the collinear region 
therefore need to be expanded at least up to this order. 

Finally, in the soft region (xi = X2 = 1), since the expansion of the distributions 
in eq. (4.4) generates additional coefficients, the function 72. is required up to 0{e^) 
where transcendental constants of weight 4 appear. In this region, we evaluated all master 
integrals by direct integration. These integrations closely resemble the integrals appearing 
in the two-loop soft function, and we used the same techniques as were applied in [56], 
expanding the resulting hypergeometric functions in e using HypExp [57]. 

The master integrals in both collinear regions and in the hard region were determined 
from their differential equations in xi and X2, using the results from the soft region as 
boundary condition. The results for all ten new master integrals can be expressed in 
terms of harmonic polylogarithms [58] and a few specific combinations of logarithms and 
polylogarithms. They are listed in the appendix. 

Insertion of the master integrals and Laurent expansion of the overall factors according 
to eq. (4.4) yields the integrated antenna functions X{xi,X2,e). 

5. Results 

The initial-initial antenna functions are listed in Tables 2-4. The full results for their 
integrated forms are too lengthy to be quoted here, a separate FORM file containing those 
is attached with the arXiv submission of this paper. In this section, we give at most 
two terms of their Laurent expansion around e = below. From each antenna function, 
we have omitted a common factor ((/^)~^^. This factor differs from the convention used 
in [42], where a common factor (512)^^'^ had been taken out. It is however in line with the 
final-final antenna functions derived in [23] and the initial-final antenna functions derived 
in [40], which are all normalised to (g^)"^*^. 

The leading pole structure of those antenna functions containing a singularity from two 
soft gluons or from a soft quark-antiquark pair can be predicted from infrared factorisation, 
and is universal between the final-final, initial-final and initial-initial configurations (in the 
former two cases, a single antenna function sometimes contains multiple soft gluon limits, 
which are then additive). For two colour- unconnected (abclian) soft gluons, we expect a 
leading pole of for two colour-connected (non-abelian) soft gluons a leading pole of 
3/4/e^ and for a soft quark-antiquark pair a soft pole of— 1/12/e^. A similar pattern can be 
seen in the integrated one-loop antenna functions [23,40,41], where the difference between 
abelian and non-abelian gluonic contributions can be explained by the soft-gluon current 
at the one-loop order [9] . The one-loop soft-gluon current contribution cancels between the 
double real radiation antennae and the one-loop antennae [23,24]. 

5.1 Quark-antiquark antenna functions 

The leading poles of the integrated initial-initial antenna functions arising from crossings 
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of the quark-antiquark antennae are as follows: 



Ai2 = Al{q,g,g,q) 



1 



^S{1 - xi)6il - X2) + ^<5(1 - X2)(l + XI - 2I?o(a;i)) 



+ -<5(l-Xi)(l+X2-2Po(x2)) 



Ai3 = Al{q,g,g,q) 

= ^4 (1 - 2x2 + 2xi) (5(1 - xi) + O (e-2) , 
Ai4 = Al{q,g,g,q) 

= ii (1 - 2x2 + 2xi) 5(1 - xi) + O (e-2) , 

= (1 - 2x1 + 2x?) (1 - 2x2 + 2xi) + O (e-i) , 

Av2 = Al{q,g,g,q) 

= ^6{1 - xi)S{l - X2) 

+4 [S{1 - X2) (1 + XI - 2Po(aJi)) + (^(1 - xi) (1 + X2 - 2Po(x2))] 
= (1 - 2x2 + 2xi) 6{1 - xi) + O (e-2) , 

^34 = ^[4(9,5,5,5) 

= 11 (1 - 2x1 + 2x?) (1 - 2x2 + 2x^) + O (e-^) , 
Bi2=Bl{q,q',q',^) 

-^5(l-XiM(l-X2) 

--^6(1 - xi)5(l - X2) - ^6{1 - X2) (1 + XI - 2Vo{xi)) 



1 



- — (5(l-Xi)(l+X2-2Po(x2)) 



Bis=B'i{q,q',q',q) 



1 



5(1 -xi) -(X2 + I)i/(0;X2) + 



(I-X2) (4x| + 7x2 + 4)\ 



24x2 



) 



+0 e 



fi34 = e2(g, 



Ci2 = C2{q,q,q,'q) 



(5.10) 



48 (1 - X2) 



((1 + xl) [6H{0, 0; X2) + 6H{1, 0; ^2) + tt^] 
+ 3(5 - 2xl)H{0; X2) + 3(1 - X2){8 - Txs)) 



+0 {e') , 
Ci3 = Cl[q,q,q,q) 



(5.11) 



<^(l-:ri) 
24(1 + X2) 



((1 + xl) [l2iJ(-l, 0; X2) - 6H{Q, 0; X2) + tt^] 



12(1 -xi)- 6(1 + a;2)^if(0;x2)) 



C23 = C2{q,q,q,q) 

C24 = C2{q,q,q,q) 
= o (e°) . 



+ O(e0), (5.12) 



(5.13) 



(5.14) 



We observe that B34, C23 and C34 are finite, as expected. Moreover, the leading pole 
structure of A12 (non-abclian), A12 (abelian) and B12 is as predicted from the double soft 
factorisation. B12, B13 and S34 were computed previously in [42]. 

5.2 Quark-gluon antenna functions 

The leading poles of the integrated initial-initial antenna functions obtained from crossing 
quark-gluon antennae read: 

= ^^S{l-xi)6{l-X2) 



1 



^S{1 - xi)S{l - X2) + ^6{1 - X2) {1 + xi- 2Vo{xi)) 



1 



+-6{1 -xi)[xi-X2 + 2 - Voix2) 



X2 



+ 0(e-2) 



(5.15) 



T^i3 = T^4{q,9,9,9) 



-jS{l-Xi)d{l-X2) 



+- 



(5(l-X2)(l + xi-2I»o(iCi)) 



+ 2(5(1 -xi)(xi-X2 + 2 - Vo{x2) 

X2 



+ O (6-2) , 



(5.16) 



V22, = 'Dl{q,g,g,g) 
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J (1 - 2x2 + 2xi) (5(1 -xi) + l{l- 2x1 + 2x?) (5(1 - X2) 



+ (6-2) , 

^^24 = 'l^'i{q,g,g,g) 



^ (1 - 2x2 + 2x1) - + ^ (1 - 2x1 + 2x^) 5{1 - X2) 



+0 (6-2) , 

1 (2 - 2x2 + xl) (5(1 - xi) 



Eli = £l{q,q' ,q' ,g) 



2X2 



+ O (6-2) , 



+ - 



^24 = £l{q,q',q\g) 



--(5(1-Xl)5(l-X2) 

-^(5(1 - xi)5(l - X2) - ^(5(1 - X2) (1 + XI - 2Po(xi)) 

, _ . f 6x2H{(}; X2)(l + X2) - 8x1 + a^j - 5x2 + 8 
+ '^^ 24x2 

^\o(6-), 

+ 0(6-i), 

(1 - Xl) (4xf + 7x1 + 4) 



+ q'Do{x2) 

£23 = £l{q,q',q',g) 

1 X1X2 (1 + xiX2)^ (xi + X2 — 2) 



(xi + X2) 



(5(1-X2) I -(l + xi)i/(0;xi) + 



+ 



(2 - 2x1 + xf) (1 - 2x2 + 2x2) 



4x1 



24x1 
+ 0(6-^), 



£12 = £l{q,q',q,g) 

1 (2 - 2x2 + xj) ,5(1 - Xl) 

6^ 4X2 



+ O (6-2) , 



5i4 = £4{q,q',q',g) 
^23 = £l{q,q',q',g) 

_ 1 X1X2 (1 + XiX2)^ {x\ + x| - 2) 



(xi + X2) 



+ 0(6-l), 



^24 = £2{q,q',q',g) 



(5.17) 



(5.18) 
(5.19) 



(5.20) 



(5.21) 



(5.22) 



(5.23) 



(5.24) 
(5.25) 
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(2xf + X2x\ - + 2x\x\ - X2x\ - 8x|xf + 2x2Xx + 8x1) 



(G(-X2; Xl) - xi)) + {2x\ + xsxi - l) x?il(0; X2) 



+ 



(1-^2) 



(I2x? + 36x2X? - 6x? + 41x^x^ - 19x2xJ 



6(xi +X2)^xf 

- x^ + IQx^xf - 29xix^' - 5x2X*i' + \x\x\ - 20x^xf - 26xixf 
+ 15x2xf - <ox\x\ - mx\x\ + Tix\x\ + 10x2xf - 2xf - 12x^xf 
+ 86xix? + 50x^x? - 34x2xf + 32x^x? + 62x^x? - 118x 



2„,2 
2-^\ 



+ 24x^X1 - 132xixi - 48x|) 



+0 (e") 



(5.26) 



Again, the leading pole structure of P12 (non-abclian), X>i3 (abelian) and £^14 is as expected. 
£"12, ^14, ^23 5 ^24 were computed already in [42]. 

5.3 Gluon-gluon antenna functions 

Finally, the leading poles of the integrated initial-initial gluon-gluon antenna functions are: 



1 3 

?4 
1 



--5il-xi)6{l-X2) 



^6(1 - xi)(5(l - X2) + ^S{1 - X2) (^x? - XI - ^ + 2 - Vo{xi) 

3/1 

+ -5(1 -Xl) xi-X2 +2-Vo{x2) 

Z \ X2 



+ 0(e-2) 



-^13 = ^{9,9,9,9) 



-^5{l - x{)5{l - X2) 



1 



-Xl) ( X2 - X2 - — 2 - Po(a::2) J 

\ X2 J 



25{l 



+2(5(1 - X2) Xl - Xl 



Xl 



+ 2-Vo{xi) 



+ O (6-2) , 



Q12 = 04(9, q,q, 9) 



+- 



--5(1-Xi)(5(l-X2) 
^5(1-Xi)(5(l-X2) 

, f 6x2-H"(0; X2)(l + X2) - 8x1 + - 5x2 + 8 

^'^V 24x2 



+ g2^0(^2) 



(5.27) 



(5.28) 
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. ,^ , / 6xiif(0; xi){l + Xl) - 8xf + xf - 5xi + 8 
+ - i 



Qi3 = Ql{g,q,q,g) 

1 (2-2x2 + xi)(5(l-xi) 
e3 2x2 

1 (2-2x2 + xi)(^(l-xi) 



2X2 



+ 0(e-2) 



+ O (6-2) , 



^34 = ^4 (5> ^> 9) 



12(1 + X1X2) (xgxf + X2X1 — ix^xl + xf + X2) 



G12 = Ql{g,q,q,g) 



(xi + X2) 



+0 (.-■), 



*(1 - X.) (i (1 + .0 ;^(0; + (l-x.)K + 7x.+4) j 



+0 (.-'), 

Q12. = G4{9,q,Q,9) 

1 (2-2x2 + xi)(5(l-xi) 
4x2 



+ 0(e-2) 



^34 = 02(9, q,Q, 9) 

12(1 + X1X2) (x2x| + xfxf — 4x2xf + xf + X2) 



(xi + X2) 



'Hi2=nl{q,%q',q') 



(1 + X1X2) (x|x| + Xgxf - 4x|xf + xf + x|) 



3 (xi + X2) 
nn = 'Hl{q,q,q',q') 

1 (2 - 2x1 + xj) (2 - 2x2 + xl) 

Ax 1X2 



+ O (eO) 



+ 0(e-). 



(5.29) 
(5.30) 
(5.31) 
(5.32) 



(5.33) 
(5.34) 

(5.35) 
(5.36) 
(5.37) 



The leading pole structure of J^i2 (non-abclian), J^i3 (abelian) and ^12 is as predicted from 
the double soft factorisation. ?ii2 and l-Li^ were computed already in [42]. 



6. Conclusions 

In this paper, we have completed the integration of double real radiation antenna functions 
in the kincmatical situation where both radiator partons are in the initial state. A subset 
of our results, restricted to antenna functions containing a secondary fermion pair, was 
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obtained earlier in [42]. With the results presented in this paper, combined with the final- 
final [23], initial- final [40] and virtual one-loop initial-initial antennae [41], all antenna 
functions required for the calculation of jet cross sections at hadron colliders are now 
available in unintegrated and integrated form. The implementation of the corresponding 
antenna subtraction terms into a numerical parton-level event generator program is thus 
readily feasible. 

Using the general method completed in this paper, NNLO QCD calculations for ex- 
clusive jet observables at hadron colliders are now in reach. The principal requirement 
for these calculations is the availability of the relevant two-loop matrix elements. Poten- 
tial applications, where two-loop matrix elements have already been derived, are: two-jet 
production [59], vector-boson-plus-jct production [60] and Higgs-boson-plus-jet produc- 
tion [61]. Important steps have already been taken in the calculation of two-jet production 
at NNLO, where the unintegrated subtraction terms for purely gluonic processes have al- 
ready been derived and tested for the double real radiation at tree-level [38] and the single 
real radiation at one-loop [39] , and in the incorporation of processes with quarks into this 
calculation [62]. Finally, we notice that the integrated initial-initial antenna functions Tvi 
and J^i3 derived here appear in the integrated subtraction terms of the all-gluon contribu- 
tion to two-jet production. They contribute to the cancellation of all infrared poles when 
combined with all other integrated subtraction terms and with the virtual two-loop matrix 
elements [63]. This pole cancellation has been verified [63] and can be regarded as a strong 
check on the correctness of our results presented here for these integrated antennae Tvi 
and T\2,- 
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A. Master integrals 

The phase space integrals appearing in the integration of the initial-initial four-parton 
antenna functions can be reduced to a set of 20 independent master integrals, as described 
in Section 4. Of these integrals, the following 10 were derived previously in [42]: 



l2 
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^3 
h 
h 
h 

Is 

/l4 
h5 



J [dk, 
J [dk, 
J [dk, 
I [dk, 
I [dk, 
J [dk, 
J [dk, 



The remaining 10 integrals are: 



A3 
ha 
In 
I22 
I23 
I24 
I25 
I27 
I28 
I29 



j[dk, 
j[dk, 
j[dh 
J[dk, 
J[dk, 

j[dh 
J[dk, 

J[dk 
J[dh 
J[dk, 



dh] <5(Ci)<5(C2) 
dkk] S{C\)6{C2) 
dkk] <5(Ci)(5(C2) 
dkk] S{Ci)6{C2) 
dkk] S{C\)5{C2) 
dkk] <5(Ci)^(C2) 
dkk] S{C\)d{C2) 
dkk] <5(Ci)<5(C2) 



-{kj ■Pi)ikj -ps) 



D 



kl2 



{kj -pi) 

Dkl2 

-{kj - Pa) 



D 



kl2 



Dkl2 
1 

Djk2 ' 
1 

1 



Djl2Djk2 
1 

DjlDk2 



[dkk] S{C,)6{C2) 
[dkk] SiCi)6{C2) 
[dkk]diCi)5{C2) 
[dkk]S{Ci)SiC2) 
[dkk] 5{Ci)d{C2) 
[dkk]5iCi)5{C2) 
[dkk] SiCi)SiC2) 
[dkk] S{Ci)d{C2) 
[dkk]6iCi)6{C2) 
[dkk] S{Ci)SiC2) 



1 



Dk2Djl2 

1 

Djk Dji Dji2 ' 

1 

Djl Dk2 Djk2 ' 

1 

Djl Dk2 Djl2 ' 

1 

Dk2 Djk2 Djl2 ' 
1 

Djl Djk2 Dji2 ' 
1 

Dj2 Dk2 Dkl2 ' 

1 

Djl Dj2 Djk2 ' 
1 

Djl Dj2 Dk2 ' 
1 

Dj2 Djk Dki 



(A.l) 



(A.2) 



Note that, the numeration of the integrals is not consecutive, since masters related by 
exchange of xi and X2 are listed only for one crossing in both sets. 

In the new set, the integrals /13, 122^ I23 and I25 appear always in particular combina- 
tions, which can be made explicit by introducing new master integrals M22, M23 and M25 
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by: 



722 = M22 + 
I23 = M23 + 



-/l3, 



2xiX2Sl2 

1 — 2x1 ~ 2xiX2 
2X2(1 - xl){xi + X2)S12 



13 



I25 = M25 + J' or^h3 ■ 

2X2(1 - Xf)si2 



(A.3) 



After expressing the integrated antenna functions in this new basis it turns out that the 
remaining coefficients of /13 are all proportional to e, such that /13 is required only to lower 
order in e than the other master integrals. 

A.l Hard region 

In the hard region, the new master integrals are: 



I 



13 



(^12) 



-2-2e 



Sr 



(l-Xi)-^^(l-X2)-^n 

Xl e 

2e/i ^ \-2e 



gi3ixi,X2) + 0{e^ 



(A.4) 



^16 = (S12) 

1 1 



3-2.5^(1-^1) (1-^2) 



+ - 



xi(l -x|) 

4 + - ( - 3G(-X2; Xl) + G(-l; xi) + 2H{-1; X2) - 3H{0- X2)) 

17^2 



12 



2G(l/x2; xi)i7(0; xa) - 2G(l/x2, 0; xi) - 6G(-X2, -X2; xi) 



+3G(-X2; xi)i7(-l; X2) - 6G{-X2; xi)H{0; X2) + 3G(-X2, -1; xi) 
+3G(-1,-X2;xi) +G(-l;xi)log2-2G(-l;xi)F(-l;x2) + 3G(-l;xi)ii'(0;x2) 
-2G(-1, -1; Xl) + 2G(0; xi)i7(0; X2) + 4G(0, 0; xi) - 3G(1, -X2; xi) 
-G(l; Xl) log 2 + 3G(1; xi)i7(-l; X2) - 3G(1; xi)i7(0; X2) + G(l, -1; xi) 
-2G(1, 0; Xl) - 4H{-1, -1; X2) + 6i7(-l, 0; X2) + 3H{0, -1; X2) 



-4i7(0,0;x2) 



+ 0{e), 



(A.5) 



Ii7 = (312)-^-^' Sr 



1 



(1-Xi)-^^(1-X2) 
XlX2(l — Xl) 



-2e 



1 

2? 



+- (log 2 - G(-l; Xl) + G(0; xi) + H{0; X2)) 

7^2 

_21og2 2 - — + G(l/x2;xi)F(-l;x2) - G(l/x2, -l;xi) + G(l/x2,0;xi) 

-3G(-1, -X2; Xl) + 2G(-1; Xl) log 2 - G(-l; xi)F(-l; X2) - G(-l; xi)/7(0; X2) 

+2G(-1, -1; Xl) + G(-l, 0; xi) - 2G(0; xi) log 2 - 2G(0; xi)H{0; X2) 

+2G(0, -1; Xl) - 2G(0, 0; xi) - 2H{0; X2) log 2 - 2i7(0, 0; X2) - 2i7(l; X2) log 2 



+2i7(l,-l;x2)-4i7(l,0;x2) 
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(1-Xi)-^^(1-X2) 



-2e 



-(-log2 + iJ(-l;x2)) 



TT 



+2 log^ 2 + — + G{l/x2;xi)H{-l; X2) - G{l/x2, -1; xi) + G(l/a;2, 0; xi) 
+2G(-1; Xl) log 2 - 2G(-1; X2) + 2G(0; xi) log 2 - 2G(0; xi)iJ(-l; X2) 

-3G(1, -X2; Xl) - 4G(1; Xl) log 2 + 3G(1; xi)iJ(-l; X2) - 3G(1; xi)iJ(0; X2) 
+4G(1, -1; Xl) - G(l, 0; Xl) - 2H{-1, -1; X2) + 2H{-1, 0; X2) + 2/f(0; X2) log 2 



-2//(0, -1; X2) + 2H{1- X2) log 2 - 2H{1, -1; X2) 



M22^(.i2)--^r^-)""(^-^^)"" 



+ 0(e), 



(A.6) 



1 1 /I 



xfx2 

^G(-X2;xi) - ^G(-l;xi) - G(0;xi) + ^//(0;x2) 



- ^G(-X2; xi)iJ(-l; X2) + G(-X2; xi)ii"(0; X2) + ^G(-X2, -1; xi) 
-G(-X2,0;xi) - ^G(-1,-X2;xi) - ^G(-l; xi) log2 + G(-l; xi)i^(-l; X2) 
-^G(-l; xi)i7(0; X2) + G(-l, -1; xi) - G(0, -X2; xi) - G(0; xi)H{0; X2) 
+G(0, -1; Xl) + 2G(0, 0; xi) - ^G(l, -X2; xi) + ^G(l; xi) log 2 
+^G(l;xi)i7(-l;x2) - ^G(l; xi)iJ(0; X2) - ^G(l,-l;xi) - G(l,0;xi) 
-2Hi-l; X2) log 2 + Hi-1, -1; X2) - Hi-1,0; X2) - ^i/(0, -1; X2) + HiO, 0; X2) 



23 



-2iJ(l; X2) log 2 + 2iJ(l, -1; X2) - 3H{1, 0; X2) - 522(xi, X2) 

(l-Xi)-2^(l-X2)-2^ 



+ 0(e), 



(A.7) 



-3-2e 



XlX2(l - xf)(xi + X2) 



^ Qg(-X2; Xl) - ^G(-l; xi) + ^H{0; X2) 

297r^ 1 
— G(-X2, -X2; Xl) - -G(-X2; xi)i7(-l; X2) - 2G(-X2; xi)i?(0; X2) 

3 5 1 

+-G(-X2, -1; Xl) - G(-X2, 0; Xl) - -^G{-\, -X2; xi) + Tpi.-^'^ log 2 

-^G(-l; xi)//(0; X2) + 2G(-1, -1; xi) + G(-l, 0; xi) - G(0, -X2; xi) 

-3G(0;xi)i7(0;x2) + G(0,-l;xi) - ^G(l,-X2;xi) - ^G(l;xi)log2 

+^G(l;xi)ii'(-l;x2) + ^G(l;xi)i/(0;x2) + ^G(l,-l;xi)-2i7(-l;x2)log2 
X2)H{1; Xl) + H{-1, -1; X2) - 3i?(-l, 0; xs) - 2H{0; X2) log 2 
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+2H{0; X2)H(l- xi) + -H{0, -1; X2) - 6H{0, 0; X2) + H{1; xi) log 2 
-3H{1; X2) log 2 + 3H{1, -1; X2) - 6H{1, 0; X2) + 4g23b{x2) - 2g22{xi, X2) 



+g23{xi,X2) +g25ixi,X2) 



+ 0(e), 



(A.8) 



^24 = (S12) 



-3-2e 



Sr 



2;ia;2(a;i - X2) 



-2e 



■{-G{-l;xi) + H{-l-X2)) 



+— + 2G(l/x2; xi)//(-l; X2) - 2G(l/x2, -1; xi) + 2G(l/x2, 0; xi) 
b 

-3G(-l,-X2;xi) + G(-l;xi)log2-3G(-l;xi)//(0;x2) + 2G(-l,-l;xi) 
-2G(0; xi)H{-l:x2) + 2G(0, -1; .xi) - 3G(1, -X2; xi) - G(l; xi) log 2 
+3G(l;xi)i/(-l;x2) -3G(l;xi)i/(0;x2) + G(l,-l;xi) -2G(l,0;xi) 
-3i7(-l; X2) log 2 - 2if(-l, -1; X2) + 3H{-1, 0; X2) - 3H{1; X2) log 2 



M25 = (S12) 



+3iJ(l, -1; X2) - 3i7(l, 0; X2) + 2^24(3^1; X2) 

(l-Xi)-2^(l-X2)-2^ 



+ 0(6), 



(A.9) 



-3-2e 



Sr 



X2(l - xj] 



^ Qg(-X2; xi) - ^G(-l; xi) + ^H{0; X2) 

SStt^ 1 1 
— 2F ~ 2'^^-^'^'^^^)^^-'^'^ ^2) + G(-X2; xi)if(0; X2) + -G(-X2, -1; xi) 

-G(-X2, 0; xi) - \g{-1, -X2; xi) + \g{-1- xi) log 2 - \g{-1- xi)H{0; X2) 
+2G(-1, -1; xi) + G(-l, 0; xi) + G(0, -X2; xi) - G(0; xi)H{0; X2) - G(0, -1; xi) 
-^G(l, -X2; xi) - ^G(l; xi) log 2 + ^G(l; xi)iJ(-l; X2) + ^G(l; xi)H{0; X2) 

+^G(1, -1; xi) - H{-1; X2) log 2 - H{-1; X2)H{1; xi) - H{-1, 0; X2) 

-2/7(0; X2) log 2 + 2if(0; X2)i?(l; xi) + ^H{0, -1; X2) - 5H{0, 0; X2) 



xi) log 2 - 2iJ(l; X2) log 2 + 2i7(l, -1; X2) - 5H{1, 0; X2) 



I27 = (S12) 



+45'236(a;2) -5'22(a;i,X2) + 5'23(a;i,a;2) + 525(3:1, X2) 

(l-Xi)-2^(l-X2)-2^ 



(A.IO) 



-3-2e 



Sr 



X2(l - xf) 



1 1 



^ + - ( - 2G(-X2;xi) +2G(-l;xi) - 2i/(0;x2)) 



+ 



e 

37r2 



2G(l/x2; xi)i7(-l; X2) + 2G(l/x2, -1; xi) - 2G(l/x2, 0; xi) 
+4G(-1, -X2; xi) + 4G(-1; xi)i7(0; X2) - 4G(-1, -1; xi) + 2G(0; xi)H{-l; X2) 
-2G(0, -1; xi) + 2G(0, 0; xi) + 2H{-1- X2) log 2 + 2i7(0, 0; X2) 



- 21 - 



^28 = (Sl2) 



+2H{1; X2) log 2 - 2iJ(l, -1; xs) + 6iJ(l, 0; X2) 

(l-Xi)-2^(l-X2)-2^ 



+ 0(e), 



(A.11) 



-3-2e 



+ \ (21og2 - G(-l;xi) + G(0;xi) - H{-\;x2) + if(0;x2)) 



27r2 

-41og^ 2 - — - 3G(-1, -X2;xi) + G(-l; xi)//(-l; X2) - G(-l; xi)i7(0; X2) 
+2G(-1, -1; xi) + 0; xi) - 4G(0; xi) log 2 + 2G(0; xi)if(-l; X2) 

-2G(0; xi)i7(0; X2) + 2G(0, -1; xi) - 2G(0, 0; xi) + 3G(1, -X2; xi) 
+4G(1; xi) log 2 - 3G(1: xi)i/(-l; X2) + 3G(1; xi)i7(0; X2) - 4G(1, -1; xi) 
+G(1, 0; xi) + 2i?(-l, -1; X2) - 2//(-l, 0; X2) - 4//(0; X2) log 2 + 2^/(0, -1; X2) 

-2iJ(0, 0; X2) - 4/f(l; X2) log 2 + 4/f(l, -1; X2) - 4/f(l, 0; X2) 

(l-Xi)-2^(l-X2)-2^ 



+ (^12) 



-3-2e 



5r 



X2(l - xf) 



- + \ (-2G(-X2;xi) + G(-l;xi) + G(0;xi) + if(-l;x2) - i/(0;x2)) 

W + G(-l, -X2; xi) + 4G(-1; xi) log 2 - 3G(-1; xi)i7(-l; X2) 
+3G(-1; xi)if(0; X2) - 2G(-1, -1; xi) + G(-l, 0; xi) + 2G(0, -X2; xi) 
-2G(0, -1; xi) - 3G(1, -X2; xi) - 4G(1; xi) log 2 + 3G(1; xi)i7(-l; X2) 
-3G(1; xi)i7(0; X2) + 4G(1, -1; xi) - G(l, 0; xi) + 2i7(-l; X2) log 2 
-2iJ(-l, -1; X2) + 2iJ(-l, 0; X2) + 2iJ(l; X2) log 2 - 2/^(1, -1; X2) 



+2i7(l,0;x2) 

-^29 = (•812)"^ 



(1-Xi)-^^(1-X2)-^^ 

xiX2(l + a;i)(l + X2) 



(A.12) 



+ - (-21og 2 + G(-l; xi) - G(0; xi) + X2) - if(0; X2)) 

27r2 

+41og2 2 + + 3G(-1, -X2; xi) - G(-l; xi)i?(-l; X2) + G(-l; xi)il(0; X2) 
-2G(-1, -1; xi) - G(-l, 0; xi) + 4G(0; xi) log 2 - 2G(0; x{)ii{-\; x^) 
+2G(0; xi)i?(0; X2) - 2G(0, -1; xi) + 2G(0, 0; xi) - 3G(1, -X2; xi) 
-4G(1; xi) log 2 + 3G(1; xi)if(-l; X2) - 3G(1; xi)iJ(0; X2) + 4G(1, -1; xi) 
-G(l, 0; xi) - 2iJ(-l, -1; X2) + 2iJ(-l, 0; X2) + 4iJ(0; X2) log 2 - 2/^(0, -1; X2) 



+2i/(0, 0; X2) + 4/7(1; X2) log 2 - 4/7(1, -1; X2) + 4i/(l, 0; X2) 



- 22 - 



x,x,{l - xl){l - xl) 

^ + ^ (- log 2 + 3G{-X2;xi) - G(-l; xi) - 2G(0; an) - H{-1- X2) + H{0; X2)) 

+41og2 2 - — + log 2 (-4G(-1; xi) + 4G(0; xi) - 4i7(-l; X2) + 4i/(0; X2)) 
+12G{-X2, -X2;xi) - 6G{-X2;xi)H{~l; X2) + 6G{-X2;xi)H{0; X2) 
-6G{-X2, -hxi) - 6G{-X2,0;xi) - 6Gi-l, -X2;xi) + 6G{-l;xi)H{-l;x2) 
-6G(-l;xi)F(0;x2) + 4G(-l,-l;xi) + 2G(-l,0;xi) -6G(0, -X2;xi) 
+2G(0, -1; xi) + 4G(0, 0; xi) + 4i?(-l, -1; X2) - 4H{-1, 0; X2) - 4H{0, -1; X2) 



+4H{0,0-X2) 



+ 0(6) . 



(A.13) 



In the above equations, we introduced the functions: 
z{xi,X2) = -^/l + 4x1X2 + 4a;| , 



5i3(a;i,X2) 
522(3:1, X2) 



1 



z{xi,X2) 



, 2;(xi,X2) + 1 - 2X2\ . fz{xi,X2) + l 
log 1 ^ TTTT— I + log 



z{xi,X2) — 1 + 2X2, 



2;(xi,X2) - ly. 



(A.14) 
(A.15) 

(A.16) 



= — + -log'^X2 + log2 log 



12 



+log 



1 - X2 
X2 



^-^] - 21ogX2 log(l + X2) + Jlog(l + X2)^ 
1 + X2 / 2 

log(z - 1) - ^log^ ( J+t) I0g(x2)l0g(l - 2X2 + Z) 



-log(l - X2)log(2: - 1 + 2x2) + log ( ^y^") log I !^2X2 



+Li2 

X2 



4X2 



2-2x2 



-Li2 
-Li2 



1-z 



+ U2 



2 - 2x2 

1-2x2 + 2 
2X2 



1 - z 

2X2 



+ Li2 



+ Li2 



2X2 



g23{xi,X2) = dyi {giz{yi,X2) - 513(1, X2)) 



2-2x2 
1-2x2-2 
2-2x2 



2X2 



log 2 log (l + + log (2 + ^ j log(2(l + X2)) 

+ log(2x2) (log(l + X2) - log(l + 2x2)) - log(4x2) log(l + 2x2) 

'1 + 2x2 + 2^ 



(A.17) 
(A.18) 



+ log(-l + 2x2 + 2) log 

+ log(l + z) (log 2 + log(x2) - log(l + 2X2 + z)) 

+ log(l - 2X2 + Z) (log 4 + log(x2) - log(l + 2X2 + z)) 



-23- 



+ log ( — TT^ ) log 



2{z-l)x^J °\l + 2x2 + z 



1+X2 

5236(a;2) = (i?(-l; y2) - log 2 - 2H{0; ya)) (A.20) 

TT^ 1 

= y + log(2a;2) log(l + 2x2) + ^U2{-1 - 2x2) + Li2(-2x2) , (A.21) 

g23cix2) = dy2 - YTm) ^^^^'■^^^ (A.22) 
= log2 + ^log32 + ^log3-21og221og3 + ^log33 + ^log(x2) 

~ log^ 21og(x2) - ^ log=^(x2) + ^ log(l + 3:2) + ^ log^ 21og(l + X2) 

1 1 

--log21og(x2)log(l + X2) - -log2(x2)log(l + X2) + 21og21og2(l + X2) 

7 1 TT^ 

+ - log(x2) log2(l + 3^2) - log^(l + X2) + Y log(l + 23^2) 

log^ 2 log(l + 2x2) + ^ log 2 log(x2) log(l + 2x2) 

1 3 
+- log2(x2) Iog(l + 2x2) - - Iog21og(l + X2) log(l + 2x2) 

-2l0g(x2)l0g(l +X2)log(l +2x2) -log^(l +X2)log(l + 2x2) 

3 7 
+- log 21og2(l + 2x2) + log(x2) log2(l + 2x2) + - log(l + X2) log2(l + 2x2) 

- log=^(l + 2x2) + ^ log 2 log (^^Ig) log(2 + 2x2) 
+\ log(l + X2) log (l±g) log(2 + 2x2) + I log(x2)Li2 (3^) 
-| log(l + X2)Li2 f — ) + Iog(x2)Li2 



2 °' ' \l + X2j " \1 + 2X2 

-log(l + 2..)Li,(^)+llog2L.,(^) 



-24- 



1 + 2x2 

2 + 2x2 



- log(l + 2x2)Li2 (3^^^) + ^ log 2 Li2 

-2Li3 - Lis - ^Li3 (\) - Li3 (^] + Li3 (-1 - 2x2) 



^ 2J V 3/ 2 V4/ V4. 

+|Li3 (-2x.) + Li3 + L.3 (^) - iLi3 (:^) 

Ji yi - X2 

log 2 1 jQg2j^-^ _|_ _ YQg2 log(l - a;2) + log(l + xi) log(xi - X2) 



2 2 

'1 + X2 



+ log(l - X2) log(l + X2) - log(xi - X2) log(l + X2) + Li2 



2 



/^l 1 — 2X2 

dyi gniyi, X2) i_^y^ (A.26) 

= + + log (-1 + ^) log(2 - 2^2) - log' (2 - 2x2) 

+ log(2 - 2x2) log(x2) + ^ log(x2) log(4x2) - log(2 - 2x2) log(z - 1) 
+ log(2x2) log{z - 1) - log(2x2) log — "'"^ 

+ log(2: - 1) log(l - 2x2 + 2;) - ^ log'(l - 2X2 + z) 
+ log(2 - 2x2) log(-l + 2x2 + z)- \og{z - 1) log(-l + 2X2 + z) 
+ log(l - 2x2 + z) log(-l + 2x2 + z) - - log'(-l + 2X2 + z) 

_Li2 - Li2 (~--^] + Li2 (~-^] +uj 



X2 J \ ^ - X2 J \ 1 - X2 J V 2(1-X2) 



with z = z(xi, X2). 

A. 2 CoUinear xi region 

In the coUinear xi region, we can expand around X2 = 1, and the master integrals read: 
7i3 = 0((1-X2)) , (A.28) 



-25- 



he = (su) Sr — 

1 1 



og2 



i + ^(2 log 2 - 2H{-1; Xi))+(^-2 log'' 2 + ^ + 2H{-1; xi) log 

-2H{-1, -1; xi) + 4H{0, 0; xi) - 2H{1; xi) log 2 + 2il(l, -1; xi) + 4il(l, 0; a;i)^ 

(33 4 
^Cs + 3 log' 2 + log2 2 xi) + 3i/(l; xi)) + log 2(2ii'(-l, -1; xi) 

+4iJ(-l, 0; xi) - 8H{0, 0; xi) - 2H{1, -1; xi) - 4H{1, 0; xi) + 4H{1, 1; xi)) 
+^7r2(- log 2 - xi) + H{-1; xi) - H{0; xi)) - 2H{-1, -1, -1; xi) 
-4H{-1, 0, -1; xi) + 8H{-1, 0, 0; xi) + 8iJ(-l, 1, 0; xi) + 4H{0, -1, 0; xi) 
+8H{0, 0, -1; xi) - 16iJ(0, 0, 0; xi) - 12/f(0, 1, 0; Xi) + 2H{1, -1, -1; xi) 
+4i7(l, -1, 0; xi) + 4H{1, 0, -1; xi) - 16H{1, 0, 0; xi) 

1, -l;xi) - 12H{1, l,0;xi)j + - xa)) 



(A.29) 



, . .-3-2. . (l-xi)-^-(l-X2r^- 
= ^'''^ xi(l-xi) 



~2^ + 7 (log 2 - H{-l;xi) + H{0;xi)) + ( " y " 21og' 2 + log 2(i/(-l; xi) 
-2Ji-(0; xi) - H{1; xi)) - ^/(-l, -1; xi) + H{-1, 0; xi) 
+2H{0, -1; xi) - 2H{0, 0; xi) + -1; xi) - H{1, 0; xi)^ 

+e(c3 + I log' 2 + ^ log2 2 (-//(-I; xi) + 4H{0; xi) + 3J?(1; xi)) 

+ log 2(iJ(-l, -1; xi) - H{-1, 0; xi) - 2H{0, -1; xi) + 4H{0, 0; xi) + 2iJ(0, 1; xi) 

-H{1, -1; xi) + 3H{1, 0; xi) + 2i/(l. 1: xi)) + y (F(-l; xi) + H{1; xi)) 
-i?(-l,-l,-l;xi) + -1,0; xi) + 0,-1; xi) - i/(-l, 0, 0; xi) 

+2i7(0, -l,-l;xi) -2i7(0, -l,0;xi) - 4if(0, 0, -1; xi) + 4if(0, 0, 0; xi) 
-2H{0, 1, -1; xi) + 2H{0, 1, 0;xi) + i/(l, -1, -l;xi) - H{1, -1, 0;xi) 

-3i7(l, 0, -1; xi) + 3H{1, 0, 0; xi) - 2iJ(l, 1, -1; xi) + 2iJ(l, 1, 0; xi)^ 
+0((1-X2)) +0(e2), (A.30) 



I22 = (.S12) 



-3-2e 



Sr 



(1-Xi)-^^(1-X2) 



-2e 



^ - ^H{0; xi) + r - i log2 2 - H{-1, 0; xi) + 2il(0, 0; xi) 



-26- 



-H{l;xi)\og2 + H{l,-l;xi) + H{l,0;xi)^ 

+e ( - C3 + ^ log' 2 + ^ log2 2 {2HiO; xi) + 3H{1; xi)) 

+ log 2 0; xi) + 2H{0, 1; xi) - -1; xi) + 0; Xi) + 2H{1, 1; xi)) 

+y xi) - xi)) - -1, 0; xi) - 0, -1; xi) 

+3i7(-l, 0, 0; xi) + 2H{-1, 1, 0; Xi) + 2H{0, -1, 0; Xi) - 4/f(0, 0, 0; xi) 
-2/^(0, 1, -1; xi) - 2if(0, 1, 0; xi) + H{1, -1, -1; xi) + H{1, -1, 0; xi) 
0, -1; xi) - 3H{1, 0, 0; xi) - 2H{1, 1, -1; xi) 

-2ii-(l,l,0;xi)^ +0((1-X2)) 



723 = 0{{1-X2)) , 



(A.31) 
(A.32) 



, . ^-3-2. (l-Xi)-^-^-(l-X2)-^- 
-'24 = — (S12) "-T 



Xi 



^(log 2 - H{-1; xi)) + ( - ^ log' 2 + log 2 {H{-1; xi) - 277(0; xi) - 27f(l; xi)) 
+27f (0, -1; xi) - 77(-l, -1; xi) + 277(1, -1; xi)^ 

+e Q log3 2 + i log2 2 (-iJ(-l; xi) + 677(0; xi) + 677(1; xi)) + log 2(77(-l, -1; xi) 

-27f (0, -1; xi) + 47f (0, 0; xi) + 477(0, 1; xi) - 277(1, -1; xi) + 477(1, 0; xi) 

, 27r2 

+47^(l,l;xi)) +— F(-l;xi) -77(-l,-l,-l;xi)+277(-l,0,0;xi) 
+27f(-l,l,0;xi) + 27r(0,-l,-l;xi) - 47^(0, 0, -1; xi) - 47f(0, 1, -1; ; 

+277(1, -1, -l;xi) - 47f(l,0, -l;xi) - 47f(l, 1, -l;xi)^ 



+0((1-X2)) 
I25 = 0{{1-X2)) , 

hr = (512)"'"^'' 'S'r 
1 



+ 0(e'), 

(l-xi)-^^(l-X2)-^^ 
(1 - X?) 



+ log^2 + 21og2(77(0;xi) + 77(l;xi)) + 



(A.33) 
(A.34) 

27f(0,-l;xi) + 27f(0,0;xi) 



-277(1, -1; xi) + 277(1, 0; xi)J +€[6(3 - log" 2 - Slog^ 2 (77(0; xi) + 7f (1; xi)) 

+ log 2(27f(0, -1; xi) - 67^(0, 0; xi) - 47f(0, 1; xi) + 27f(l, -1; xi) 

27r2 

-677(1, 0; xi) - 477(1, 1; xi)) log 2 - (log 2 + 77(0; xi) + 77(1; xi)) 
-277(0, -1, -1; xi) + 277(0, -1, 0; xi) + 677(0, 0, -1; xi) - 677(0, 0, 0; xi) 



-27- 



+4H{0, 1, -1; xi) - AH{0, 1, 0; xi) - 2H{1, -1, -1; xi) + 2H{1, -1, 0; xi] 
+6H{1, 0, -1; xi) - 6H{1, 0, 0; xi) + 4iJ(l, 1, -1; xi) - iH{l, 1, 0; xi) 



+0((1-X2)) 
-^28 = (Sl2)~^~^^ 'S'r 



(l-xi)-^^(l-3;2)-'^ 

Xl(l - Xl) 



-2i7(0; xi) - H{1; xi)) - -1; xi) + H{-1, 0; xi) 

+2ii-(0, -1; xi) - 2//(0, 0; xi) + -1; xi) - 0; xi) 



+e( Cs + ^ log^ 2 + ^ log2 2 xi) + 4H{0; xi) + 3iJ(l; xi)) 



+ log 2(^i7(-l, -1; xi) - //(-1, 0; xi) - 2i7(0, -1; xi) 

+4i/(0, 0; xi) + 2H{0, 1; xi) - 17(1, -1; xi) + 3H{1, 0; xi) + 2iJ(l, 1; xi)) 

+ y (i7(-l;xi) + H{l;xi)) - H{-1, -1, -l;xi) + H{-1, -l,0;xi) 
+i7(-l,0, -l;xi) - if(-l,0,0;xi) + 2i?(0, -1, -l;xi) - 2i7(0, -l,0;xi) 
-4i7(0, 0, -1; xi) + 4H{0, 0, 0; xi) - 2H{0, 1, -1; xi) + 2H{0, 1, 0; xi) 
-1, -1; xi) - -1, 0; xi) - 3if(l, 0, -1; xi) + 3if(l, 0, 0; xi) 



2H{1, 1, -1; xi) + 2H{1, 1, 0; xi)^ + - xs)) 

(l-Xi)-2^(l-X2)-2^ 



+ (S12) 



-3-2e 



5r 



1-x? 



+ log2 2 + ^ + log 2 {2H{0; xi) + 2H{1; xi)) - 2i/(0, -1; xi) 



+2if(0, 0; xi) - 2H{1, -1; xi) + 2i7(l, 0; xi] 

+e(^6C3 - log' 2 - Slog' 2(/^(0; xi) + H{1; xi)) + log 2(2iJ(0, -1; xi) 

-6i/(0,0;xi) -4i7(0, l;xi) + 2i7(l,-l;xi) -6i?(l,0;xi) - 4i7(l, 1; xi)) 



(A.35) 



TT 



+y (-2 log 2 - 2HiO- xi) - 2i?(l; xi)) - 2^/(0, -1, -1; xi) + 2J/(0, -1, 0; xi) 
+6H{0, 0, -1; xi) - 6H{0, 0, 0; xi) + 4H{0, 1, -1; xi) - 4F(0, 1, 0; Xi) 
-2H{1, -1, -1; xi) + 2H{1, -1, 0; Xi) + 6iJ(l, 0, -1; xi) - 6i7(l, 0, 0; xi) 



+4i/(l, 1, -1; xi) - 4Hil, 1,0; xi) j + - xs)) 



(A.36) 



-28- 



'29 



{S12) 



-3-2e 



St 



(l-Xi)-^-^^(l-X2) 

2x1 



-l-2e 



^ + 1 (-4 log 2 + 4iJ(-l;a;i) - AH{0;xi)) + (^21og2 2 - ^ 

+ log 2 (-4/f(-l; xi) + 4/^(0; xi)) + 4//(-l, -1; xi) - 4i/(-l, 0; xi) 

-41^(0, -1; xi) + 4/7(0, 0; xi)^ + - UCs + ^ (log 2 - i7(-l; xi) + H{0; xi)) 

+41og2 (-i?(-l, -1; xi) - H{-1, 1; xi) + H{0, -1; xi) + HiO, 1; xi)) 
+4/7(-l, -1, -1; xi) - 4iJ(-l, -1, 0; xi) + 4H{-1, 1, -1; xi) - 4iJ(-l, 1, 0; xi) 

-4i7(0, -1, -1; xi) + 4iJ(0, -1, 0; xi) - 4iJ(0, 1, -1; xi) + 4/f(0, 1, 0; xi) 



+0((1-X2)) 



(A.37) 



A. 3 CoUinear X2 region 

In the coUinear X2 region, we can expand around xi = 1, and the master integrals read: 

, . ^-3-2.c, (l-Xi)-^-(l-X2)-^^ 

^ (- log 2 + H{-1; X2) - 2H{0; xs)) 

+ Q log2 2 + log 2 {-H{-1; X2) + 47/(0; X2) + 2H{1; X2)) + ^ + ^^(-1> -1; 3^2) 
-2H(-1, 0; X2) - 4.H{0, -1; X2) + 6i7(0, 0; X2) - 2i7(l, -1; X2) + 41/(1, 0; X2) 

+e (^4(3 - ^ log^ 2 + log2 2 Qi^(-1; ^2) - 4/1(0; X2) - 3il(l; X2)^ 

+ log 2( - il(-l, -1; X2) + 2il(-l, 0; xs) + 411(0, -1; xs) - 1011(0, 0; X2) 
-411(0, 1; X2) + 211(1, -1; X2) - 811(1, 0; xs) - 41^(1, 1; X2)) 
27r2 

-— (log 2 + 11(0; X2) + 211(1; xs)) + 11(-1, -1, -1; X2) - 211(-1, -1, 0; X2) 
-211(-l,0,-l;x2) +41^(^1,0,0; X2) + 211(^-1,1, 0;x2) - 41/(0, -1, -1; X2) 
+611(0, -1, 0; X2) + 10ir(0, 0, -1; X2) - 14ir(0, 0, 0; X2) + 411(0, 1, -1; X2) 
-811(0, 1, 0; X2) - 211(1, -1, -1; X2) + 411(1, -1, 0; X2) + 811(1, 0, -1; X2) 

-1211(1, 0, 0; X2) + 411(1, 1, -1; X2) - 811(1, 1, 0; X2) 



+0((l-xi)) 



(A.38) 



'16 



{312)-'-'' Sr 



(1-Xi)-^^(1-X2) 

1-xi 



-2e 
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-4 + -(log2-iJ(-l;x2)) + 



1 



log2 2 + iJ(-l;x2)log2 + 



1 



27r2 



-H{-1, -1; X2) + 2i?(0, 0; X2) + 21^(1, 0; X2) J + 6(6(3 + ^ V 2 

log2 2/f(-l; ,X2) + log2 -1; X2) - 2H{0, 0; X2) - 2H{1, 0; X2)) 



+ 



27r2 



(- log 2 + H{-1; X2) - H{0; X2) - H{1; X2)) - H{-1, -1, -1; X2) 



+2H{-1, 0, 0; X2) + 2H{-1, 1, 0; X2) + 2H{0, -1, 0; X2) + 2/^(0, 0, -1; X2) 
-6H{0, 0, 0;x2)- 4iJ(0, 1, 0; X2) + 2/^(1, -1, 0; X2) + 2H{1, 0, -1; X2) 



fl7 



-6i/(l, 0, 0; X2) - 4i7(l, 1, 0; X2)) + - a;i)) 



(A.39) 



2x1 



+ ^/f(0; X2) + (^2 log^ 2 + 4iJ(l; X2) log 2 - y - 4H{0, 0; X2) - 4if(l, -1; X2) 

-2H{1, 0; X2)) + - 4C3 - ^ V 2 + V 2 (-4iJ(0; X2) - 8/^(1; X2)) 

+ log 2 (-8iJ(0, 1; X2) + 4H{1, -1; X2) - 4H{1, 0; X2) - 12i7(l, 1; X2)) 

+— (2 log 2 + 2/7(0; X2) + 5H{1; X2)) + 8H{0, 0, 0; X2) + 8F(0, 1, -1; X2) 
+417(0, 1, 0; X2) - 4i7(l, -1, -1; X2) + 417(1, -1, 0; X2) + 477(1, 0, -1; X2) 

+677(1, 0, 0; X2) + 1277(1, 1, -1; X2) + 477(1, 1, 0; X2) 



+0((l-xi)) 



+ 0(e2). 



(A.40) 



M22 = (si2)-'-'^5r 



(l-xi)-^^(l-X2) 

X2 



\-2e 



_L + J_(_iog2 + 77(-l;x2)) 



11 1 

+- log2 2 - -77(-l; X2) log 2 - y + -77(-l, -1; X2) - 77(0, 0; X2) - 77(1, 0; X2) 
+e ^ - 3C3 - ^ log3 2 + i log2 2 77(-l; X2) 
+ log 2 (^—Hi-l, -1; X2) + 77(0, 0; X2) + 77(1, 0; X2) 

+y (log 2 - 77(-l; X2) + 77(1; X2) + 77(0; X2)) + ^77(-l, -1, -1; X2) 
-77(-l, 0, 0; X2) - 77(-l, 1, 0; X2) - 77(0, -1, 0; X2) - 77(0, 0, -1; X2) 
+377(0, 0, 0; X2) + 277(0, 1, 0; X2) - 77(1, -1, 0; X2) - 77(1, 0, -1; X2) 



+377(1, 0, 0; X2) + 277(1, 1, 0; X2) + 0((1 - xi)) 



+ 0(e2 



(A.41) 
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'23 



(^12) 



-3-2e 



Sr 



{l-Xi)-'-'%l-X2) 
2x2(1 +X2) 



-2e 



2 47r2 

-i7(0; X2) - — - 4i7(0; X2) log 2 - 4H{-1, 0; X2) + 4i7(0, -1; X2) - 6/^(0, 0; X2) 

6 o 



-^24 = (S12) 



AH{1, 0; X2) + 45236(X2) + 0((1 - xi)) 

(l-Xi)-2^(l-X2)-l-2^ 



+ 0(6), 



(A.42) 



-3-2€ 



Sr 



X2 



-(- log 2 + iJ-(-l; X2)) + - log^ 2 - — + log 2{-H{-l; X2) + 2H{1; X2)) 



6 



13 



7 



+i/(-l,-l;x2)-2i/(l,-l;x2) J + e - ^Cs - g V 2 

+^ Iog2 2(ii-(-l; X2) - 6i7(l; X2)) + log 2(2/7(1, 0; X2) + 2i7(0, 0; X2) - 4iJ(l, 1; X2) 

+2F(1, -1; .T2) - i7(-l, -1; X2)) + y (3 log 2 - 2//(-l; X2) + i/(l; X2)) 
+i7(-l, -1, -1; X2) - 2/i'(-l, 0, 0; X2) - 2ii'(-l, 1, 0; X2) - 2F(0, -1, 0; X2) 
-2H{0, 0, -1; X2) - 2H{1, -1, -1; X2) - 2H{1, -1, 0; X2) - 2H{1, 0, -1; X2) 



+4/7(1, 1,-I;x2)) +C»((l-xi)) 



+ 0(e'), 



(A.43) 



, _ . ^-3-2. . (l-a:i)-^-^-(l-X2)-^- 

-'25 = (Si2j Or 

2X2 



-H{0; X2) - TT^ - 4i?(0; X2) log 2 + 4iJ(0, -1; X2) - 8H{0, 0; X2) 
-4i7(l, 0; X2) + i923b{x2) + 0((1 - xi)) + 0(e) , 



(A.44) 



2X2 

+ -(2 log 2 - 2i7(-l; xs)) + ( - 2 log^ 2 - 4 log 2i7(l; X2) + tt^ + 2/7(0, 0; X2) 

+4/7(1, -1; X2) + 2/7(1, 0; X2)) + e (^y (3 + ^ log' 2 + V 2(2/7(-l; X2) 
+6/7(1; X2)) + log 2(4/7(-l, 0; X2) + 4/7(-l, 1; X2) - 4/7(0, 0; X2) + 12/7(1, 1; X2)) 

+— (-6 log 2 + 4i7(-l; X2) - 3/7(0; X2) - 5/7(1; X2)) - 4/7(-l, 0, -1; X2) 
+4/7(-l, 0, 0; X2) - 4/7(-l, 1, -1; X2) + 4i7(-l, 1, 0; X2) + 4/7(0, 0, -1; X2) 
-6/7(0, 0, 0; X2) - 4/7(0, 1, 0; X2) - 6/7(1, 0, 0; X2) - 12/7(1, 1, -1; X2) 

+ 0{e^), (A.45) 



-4/7(1, 1,0; X2)) +0((l-xi)) 



- 31 - 



'28 



(^12) 



-3-2e 



Sr 



2X2 



-2e 



-| + ^ (2 log 2 - 2/f (-1; X2) + 2H{0; 0^2)) + ^ + e (^4(3 - ^ log=^ 2 

+ log2 2 (2//(-l; X2) - 4H{0; X2) - 2H{1; X2)) + log 2(4i/(-l, 0; X2) 

+AH{-1, 1; X2) + 4i7(0, -1; X2) - 8H{0, 0; X2) - 4H{0, 1; X2) + 4H{1, -1; X2) 
^ 47r2 

-4i7(l, 0; X2)) + — (- log 2 + X2) - H{0; X2)) - 4F(-1, 0, -1; X2) 

+4/i'(-l, 0, 0; X2) - 4//(-l, 1, -1; X2) + 4F(-1, 1, 0; X2) - 4F(0, -1, ~1; X2) 
+AH{0, -1, 0; X2) + 8H{0, 0, -1; X2) - 8F(0, 0, 0; X2) + 4H{0, 1, -1; X2) 
-4H{0, 1, 0; X2) - 4H{1, -1, -1; X2) + 4H{1, -1, 0; X2) + 4H{1, 0, -1; X2) 



'29 



4/f(l,0,0;x2)^ +C»((l-xi)) 

(l-Xi)-l-2^(l-X2)-l-'^ 



(A.46) 



(^12) 



-3-2^5^ 



1 



2X2 



47r^ 



- + - (-4 log 2 + 4iJ(-l; X2) - 4iJ(0; X2)) + 2 log" 2 - — 

6 6 \ o 

+ log 2 (-4i7(-l; X2) + 4H{0; X2)) + 4/f(-l, -1; X2) - 4if(-l, 0; X2) 

-m{0, -1; X2) + 4^/(0, 0; X2)) + e - MCa + ^ (log 2 - H{-1; X2) + //(O; X2)) 

+ log 2 (-4iJ(-l, -1; X2) - 4F(-1, 1; X2) + 4F(0, -1; X2) + 4iJ(0, 1; X2)) 
+4if(-l, -1, -1; X2) - 4iJ(-l, -1, 0; X2) + 4if(-l, 1, -1; X2) - 4H{-1, 1, 0; X2) 

-4H{0, -1, -1; X2) + 4iJ(0, -1, 0; X2) - 4H{0, 1, -1; X2) + 4H{0, 1, 0; X2) 



+0((l-xi)) 



(A.47) 



The combination (1 + X2)/23 — ^25 appears in this collinear limit to order 0{e), it reads: 

(l-Xi)-l-2^(l-X2)-2^ 



(1 + X2)/23 - I25 = (S12) 

r.2 



-3-2e 



Sr 



2X2 



TT 



y-4i/(-l,0;x2) + 2^/(0, 0;x2) 

+€ ^ - 3OC3 + log 2 {8H(-l, 0; X2) - 4H{0, 0; X2)) 
27r2 



+ 



(log 2 + 4F(-1; X2) - 2H{0; X2) + //(I; X2)) 



+8//(-l, -1, 0; X2) - 8if(-l, 0, -1; X2) + 12F(-1, 0, 0; X2) 
+8H{-1, 1, 0; X2) - 4F(0, -1, 0; X2) + 4F(0, 0, -1; X2) 
-WH{0, 0, 0; X2) - 4H{0, 1, 0; X2) + 8H{1, -1, 0; X2) 



-32- 



-4i7(l, 0, 0; X2) + 8523c(x2) + - xi)) 



+ 0(e2) , (A.48) 



A. 4 Soft region 

In the soft region, xi — )■ 1 and a;2 — )■ 1. The master integrals then become: 



/l3 


= {S12)-' 


/l6 


= -(^12)- 


/l7 


= -(^12)- 


-^22 


= (^12)-^ 


^23 


= -(^12)- 


hi 


= (^12)"' 




= -(^12)- 


hi 


= -(^12)- 


hs 


= -(^12)- 


hg 


= (^12)-' 



2e(l - 2e) 



-3-2.c_(l-^ir'^(-l-^2)^-'^ 



Sr 



2e2 



.3_2, „ (l-xi)-^-^-(l-X2)-^- 


2? 

.3_2_ (1-Xi)-1-2^(1-X2) 



p 


- 4C3e' 




3 


- 4C3€' 





l-2t 



2e(l - 2e) 
(1 - - X2)-^- 
2^^ 

(l-xi)-i-2^(l-X2; 



l-2e 



Sr 



e(l - 26) 



l-^^^-4C3e^-^e^ + 0(e^) 



2i^ 

-3-2. c (l-Xl)-l-2e(l_^,)-2. 



.2,^ (1-Xir-^^(1-X2) 



-l-2e 



3 

1 - y - 4C3e^ - p + 
3 27r2 



37r^ 



(A.49) 
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